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Abstract. Magnetic field amplification is needed to accelerate cosmic cays to PeV 
energies in supernova remants. Escaping cosmic rays trigger a return current in the 
plasma that drives a non-resonant hybrid instability. We run simulations in which we 
represent the escaping cosmic rays with the plasma return current, keeping the max- 
imum cosmic ray energy fixed, and evaluate its eff'ects on the upstream medium. In 
addition to magnetic field amplification, density perturbations arise that, when passing 
through the shock, further increase amplification levels downstream. As the growth rate 
of the instability is most rapid for the smaller scales, the resolution is a limiting factor 
in the amplification that can be reached with these simulations. 



1. Introduction 



Cosmic rays (CRs) are accelerated in supernova remnants (SNRs). They can be ob 
served in situ mainly in the form of relativistic electrons, that radiate synchrotron radia- 
tion from radio- to X-ray wavelengths (Reynolds & Ellison|1992 ; Koyama et al.|1995] ). 



From the radiation signature the energies for electrons and magnetic field strengths can 
be deduced (Vin k & Laming|[2003| ). Acceleration of electrons to TeV energies seems 



ubiquitous, and diffusive shock acceleration is the preferred acceleration method. Pro- 
tons are also believed to be accelerated in this process, possibly to higher energies than 
the electrons, as they do not suffer the same amount of losses. They could leave a sig- 
nature in gamma rays, which is, however, difficult to disentangle from inverse compton 
emission of the electron population. Future observations of high-energy gamma rays by 
CTA ( |Actis et al.|2011[ ), in combination with theory, could shed light on this ambiguity. 



For efficient acceleration, a magnetic field that is strong and turbulent enough 
needs to be present upstream of the shock, such that the cosmic rays are scattered back 
and forth across the shock, gaining energy at each shock crossing. The most likely 
candidate to amplify the field is the non-resonant hybrid instability that grows most 
rapidly on small scales Bell (2004). We present simulations of this instability in the up- 
stream of a SNR shock wave. It can be shown to be responsible for creating significant 
amplification upstream of the shock. In addition to amplified magnetic field it creates 
density inhomogeneities. When the shock crosses this region of density and magnetic 
field fluctuations, additional amplification occurs downstream ( |Beresnyak||2011t |Guo 
etaL|2012D 
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Figure 1. Evolution of the upstream medium of a SNR, cropped at the bottom. 
Both the density and 5^ are shown, time progresses left to right (0, 80, 160 yrs). 
The angle that is simulated is ~ 3° wide, corresponding to a physical length scale of 
~ 2.4 X 10^^ cm at the maximum radius that is plotted (5x10^^ cm). 



2. Current driven MFA 

If cosmic rays are efficiently accelerated at the shock, it can be assumed that a certain 
fraction;^ of the shock's kinetic energy is transferred to the relativistic particle popu- 
lation. If this fraction is independent on what happens at the high-energy end of the 
cosmic ray spectrum, the escaping cosmic rays will give rise to a current that decreases 
when the maximum particle energy increases. The return current that is consequently 
triggered in the plasma can give rise to a number of current driven instabilities (for an 
overview see jSchure et al.|20T2l ). 



We are investigating the effect of this current on the magnetic field growth and 
stability of the upstream plasma in the regime of the non-resonant hybrid instability 
( |Bell|2004| ). Its maximum growth rate and corresponding wave number are: 
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In order to be able to resolve this wavenumber in simulations, the size of the grid cells 
needs to be smaller than Axg < 2n/k. 

We have run a simulation in which we let a SNR evolve in a medium that initially 
has a homogeneous density, and a background magnetic field ^o, to which we add a 
turbulent magnetic field with fluctuations on rather small scales, such that they are re- 
solved by a few grid cells. We use the MHD code from the AMRVAC framework (l van| 
der Hoist & Keppens|2007| ) on a 2.5D spherical grid, with symmetry in the 3rd dimen- 



sion. Bq is in the direction of the symmetry axis. Because of resolution requirements 
we only take a small wedge, over an angle ~ 3° and radial extent of 5 x 10^^ cm. In this 
case, the smallest grid size is approximately 5 x 10^"^ cm, which is inversely linked to 
the maximum wavenumber of the instability that we can expect to grow. The current is 
representative for particles escaping beyond 10^^ eV, with 30% of the shock's kinetic 
energy going into cosmic rays. The fastest growing modes is parallel to the zeroth- 
order magnetic field, and it drives growth of the magnetic field perpendicular to the 
zeroth order current As the growth rate increases with wave number, the limited 
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resolution in this simulation means that we do not expect to resolve amplification at 
the maximum level. The general trend though can be seen in the results presented in 
Figure [T] which shows the evolution of the density and with time. Only a section of 
the simulated region is shown, such that the shock can be seen to enter the frame only 
at the last shown time step. 

The Lorentz force resulting from the jo x pushes around the plasma, thereby 
locally amplifying the magnetic field where the density is increased, while increasing 
the length scale of the structures. The maximum growth is limited because of the reso- 
lution and the small cosmic ray current: SB/Bq grows only a factor of a few in the time 
available until the shock passes. When the growth of fluctuations is faster, the density 
perturbations can be seen to disturb the shock front. 
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Figure 2. Evolution of 5^ as a function of time (top to bottom: 0, 40, 80, 120, 
160 yrs) for different magnetic field strengths: the left column has = 5yuG, the 
right column has = lOyuG. 
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Simulations on smaller scales give larger growth rates, as the growth rate goes 
up with wave number. In Figure [2] we present the evolution of the field strength for a 
simulation of just a small part of the upstream (for radius r = [4.5 - 4.55] x lO^^cm). 
To further aid the instability we have increased the current to represent escape of CRs 
beyond 10^^ eV. The left column shows snapshots of the simulation for a background 
magnetic field strength of 5 yuG, versus a magnetic field strength in the right column of 
10 fiG. For both cases SB/Bq is on the 5% level, with maximum fluctuations up to 15%. 

Growth of the instability is faster for the stronger magnetic field case, as the jo x 
Bj_ has a larger magnitude. Because of the lower magnetic field strength in the left 
column, the wave number of the fastest growing mode is significantly smaller than the 
scale of the initial perturbations, as can be seen especially in the third row. As the 
shock comes closer, the current increases. Together with the already amplified 5_l the 
growth proceeds rapidly to larger scales and amplitude. In the right column, values of 
SB /Bo as high as ~ 35 are reached. At this level, the self-consistency of the simulation 
can obviously be questioned. As the amplitude of fluctuations increase, in reality the 
cosmic rays will be confined closer to the shock, thereby decreasing the current of 
escaping cosmic rays. Additionally, the gyro radius may be smaller than the scale of 
the fluctuations, such that the assumption of a radial zeroth order current is not strictly 
valid. It remains to be investigated what happens when this feedback is included in 
simulations. 

3. Discussion 

When the turbulent magnetic field is amplified, it confines cosmic rays more effectively 
to the shock region, thus allowing more efficient acceleration. The magnetic fluctua- 
tions grow most rapidly on the small scales, and increase in scale as the non-resonant 
hybrid instability progresses. The length scale of the field fluctuations determines the 
energy to which cosmic rays are confined. When cosmic rays are better confined, the 
current will decrease, thereby in turn lowering the growth rate of the instability. 
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